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For many physical systems which can be approximated by a classical background field plus small 
(linearized) quantum fluctuations, a fundamental question concerns the correct description of the 
backreaction of the quantum fluctuations onto the dynamics of the classical background. We inves- 
tigate this problem for the example of dilute atomic/molecular Bose-Einstein condensates, for which 
the microscopic dynamical behavior is under control. It turns out that the effective-action technique 
does not yield the correct result in general and that the knowledge of the pseudo-energy-momentum 
tensor (T M „) is not sufficient to describe quantum backreaction. 

PACS numbers: 03.75.Kk, 04.62. +v, 43.35.+d, 04.60.-m 



I. INTRODUCTION 

Many phenomena in physics can be described, to a suf- 
ficient degree of accuracy, by means of the background- 
field method, wherein one considers linearized quantum 
fluctuations on top of an approximately classical back- 
ground. For a scalar gap- less (e.g., Goldstone) mode </>, 
propagating in an arbitrary background (e.g., zeroth 
sound in inhomogeneous superfluids), the low-energy ef- 
fective action can be cast into the pseudocovariant form 
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£eff = 7; {d u <t>)G» v {d v <t>) 



(1) 



with — (d/dt, V) denoting the space-time derivatives 
(we employ the summation convention), and the space- 
time dependent tensor G^ v characterizing the back- 
ground 0, 0- Hence the propagation of <ft is com- 
pletely analogous to that of a scalar field in a curved 
space-time with the effective metric g^ depending on 
the background tensor via G^ v = Iffcff l^ffcff ' wnere 
g c g = det[gf^,] (with possibly an additional dilaton fac- 
tor). This striking analogy allows us to apply many 
tools and methods developed for quantum fields in curved 
space-times |3|, and to conclude that associated phe- 
nomena like Hawking radiation yj, super-radiance etc., 
may occur. Since the equation of motion for the scalar 
mode, Dcff^ = V° ff V° ff = 0, is equivalent to covari- 
ant energy-momentum balance, y^T^ = 0, one such 
tool is the pseudo-energy-momentum tensor, cf. Q 



(2) 



The identification of effects by analogy, described above 
(including the existence of a pseudo-energy-momentum 
tensor), is of a kinematical nature. An important ques- 
tion concerns the dynamics, that is, the backreaction 
of these quantum fluctuations onto the classical back- 
ground. Extending the analogy to curved space-times a 
bit further, one is tempted to apply the effective-action 



method (see, e.g., 3] and 0), and, since the depen- 
dence of the effective action A e s on the degrees of free- 
dom of the background i] enters via the effective metric 
= g^g(n), to determine the backreaction contribution 
to their equation of motion by 



8Aes _ SA c s Sg^a _ 1 



8r) 



drj 



(3) 



The precise meaning of the expectation value of the 
pseudo-energy-momentum tensor, (T„„), is difficult to 
grasp in general, due to the nonuniqueness of the vacuum 
state in a complicated curved space-time background and 
the ultraviolet (UV) renormalization procedure. Adopt- 
ing a covariant renormalization scheme, the results for 
(Tfiu) can be classified in terms of geometrical quantities 
(cf. the trace anomaly 0). 

However, in calculating the quantum backreaction in 
that way, one is implicitly making two essential assump- 
tions: firstly, that the leading contributions to the back- 
reaction are completely determined by the effective ac- 
tion in Eq. Q , and, secondly, that deviations from the 
low-energy effective action at high energies do not af- 
fect the (renormalized) expectation value of the pseudo- 
energy-momentum tensor, (T M „). Since the effective co- 
variance in Eq. is only a low-energy property, the ap- 
plicability of a covariant renormalization scheme is not 
obvious in general. In the following, we critically examine 
the question of whether the two assumptions mentioned 
above are justified, e.g., whether {T^ v ) completely deter- 
mines the backreaction of the linearized quantum fluctu- 
ations. In oder to address this question, we shall consider 
phonon modes in a well understood superfluid, a dilute 
Bose-Einstein condensate. 

This paper is organized as follows. In Section [HJ we 
give a brief introduction to Bose-Einstein condensates, 
and introduce a particle-number-conserving ansatz for 
the field operator. In the subsequent Section ITlTl based 
on this ansatz, the backreaction to the full current will 
be calculated, yielding a different result than that ob- 
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tained by the effective action method. The differences of 
the two expressions will be investigated in more detail in 
Sec. lIVI Since the choice of the (classical) background is 
of particular importance, we will demonstrate its influ- 
ence in Sec. [V] where we consider two different choices 
for the background yielding two different expressions for 
the backreaction. Afterwards, the failure of the effective 
action technique is discussed in more detail in Sec. IVII 
The cutoff dependence of the pseudo-energy-momentum 
tensor is addressed in Sec. lVffl As a simple example, 
we consider the influence of the backreaction contribution 
to a static ID condensate in Sec. IVIIII 



II. BOSE-EINSTEIN CONDENSATES 

In the usual s-wave scattering approximation, the 
many-particle system is described by the field operator 
in the Heisenberg picture (we set ft = m = 1) 



.8 f 



> 2 



V + g¥^> I * 



(4) 



where V denotes the one-particle trapping potential, and 
the coupling constant g is related to the s-wave scattering 
length a s via g — Aira s (in three spatial dimensions). In 
the limit of many particles JV > 1, in a finite trap at 
zero temperature with almost complete condensation, the 
full field operator ^ can be represented in terms of the 
particle- number-conserving mean- field ansatz 0, Q 



(A + x + c) i/v 7 ^, 



(5) 



with the order parameter tp c — 0(\JN), the one-particle 
excitations x — O(N ), and the remaining higher-order 
corrections f = 0{1/ \f~N). The above mean-field ansatz 
can be derived in the dilute-gas limit by formally setting 

g = 0(l/N) HEHII3; we sha11 use this formal defini - 
tion of the dilute-gas limit in what follows. Insertion of 
Eq. JSJ into (0} yields the Gross-Pitaevskii equation 
for the order parameter ip c 



9 / 



iv 2 + ^ + .g|^ c | 2 + 2.g(x t x)^c 

+g(x 2 )r c , (6) 

the Bogoliubov-de Gennes equations ^3] f° r X 

= ("5 v " + v + 2 - 9l ^ c|2 ) * + 9 ^ ' (7) 

and for the remaining higher-order corrections 
C = 0{l/y/N), neglecting 0(1/N) terms: 
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-iv 2 + y + 2. g |^ c | 2 )c + # c 2 C t 



+2ff(x f X " (x f X»^c + 5 (X 2 " (x 2 »V-c* • (8) 



Note that the additional terms 2g{x X) ano - dix 2 ) m the 
Gross-Pitaevskii equation © ensure that the expecta- 
tion value of £ = 0(l/y/~N) vanishes in leading order, 
(C) = 0(1/N), see also 0,113 Without these additional 
terms, the mean-field expansion (JSJ would still be valid 
with C = 0(1/VN), but without (C) = <D(1/N). 



III. BACKREACTION 

The observation that the Gross-Pitaevskii equation © 
yields an equation correct to leading order 0(y/~N), using 
either \ip c \ 2 or IV'cl 2 + 2 (y X.) m the first line of ©, hints 
at the fact that quantum backreaction effects correspond 
to next-to-leading order terms and cannot be derived ab 
initio in the above manner without additional assump- 
tions. Therefore, we shall employ an alternative method: 
In terms of the exact density and current given by 



- / At 



,j = - <*tv*-H.c. 



(9) 



the time-evolution is governed by the equation of conti- 
nuity for g and an Euler type equation for the current j. 
After insertion of Eq. Q , we find that the equation of 
continuity is not modified by the quantum fluctuations 
but satisfied exactly (i.e., to all orders in 1/N or ft) 



(10) 



in accordance with the Noether theorem and the U(l) in- 
variance of the Hamiltonian, cf. pj. However, if we insert 
the mean-field expansion JSJ and write the full density as 
a sum of condensed and noncondensed parts 



g = Q c + (x i x) + o(i/VN), 



(11) 



with g c — \ip c \ 2 , we find that neither part is con- 
served separately in general. Note that this split re- 
quires (Q = 0(1/N), i.e., the modifications to the Gross- 
Pitaevskii equation @ discussed above. Similarly, we 
may split up the full current [with g c v c = 3(i/i*Vti) c )] 

3 = Qcv c + ~~ (x f Vx - H.c.) + 0(1/ VN) , (12) 

and introduce an average velocity v via j = gv. This 
enables us to unambiguously define the quantum backre- 
action Q as the following additional contribution in an 
equation of motion for j analogous to the Euler equation: 



^- t j = f(J,e) + Q + o(i/VN), 

where (j = gv) 

fti, 0) = -v [V • (gv)] - g{v ■ V)v 
'lV 2 ^g 



(13) 



-f?V 



2 Ve 



V cxt - gg (14) 
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denotes the usual classical force density term. Formula- 
tion in terms of the conventional Euler equation yields 



+ v ■ V ) v = V 



V 



99- Y 



iv 2 ^ 



+q/q + o(i/Vnz). 



(15) 



The quantum backreaction Q can now be calculated by 
comparing the two equations above and expressing dj / dt 
in terms of the field operators via Eqs. I0J and © 



fL- = I/^t V 3 ^-(V 2 * t )V# + H.c. 
dt J 4 \ y ' 



1 



$ T f)VV r -^-V( fi 2 (* I ) 2 f 



(16) 



After insertion of the mean- field expansion (JSJ , we obtain 
the leading contributions in the Thomas-Fermi limit 



Q = V • (v ® j, 



v — g x v (g> 



-- V-((Vx f )®Vx + H.c.) , 



(17) 



with p x = (jfix) an( l Jv = ^(x^x)- Under the as- 
sumption that the relevant length scales A for varia- 
tions of, e.g., g and g, are much larger than the heal- 
ing length £ = (<?£>c) we have neglected terms con- 
taining quantum pressure contributions ~V 2 g and [Vg] 2 
(Thomas- Fermi or local-density approximation). These 
contributions would, in particular, spoil the effective (lo- 
cal) geometry in Eq. (JIJ . 



IV. COMPARISON WITH EFFECTIVE-ACTION 
TECHNIQUE 

In order to compare the expression (|17fl with Eq. 
we have to identify cf> and . Phonon modes with wave- 
length A ^> £ are described by the action in Eq. Q in 
terms of the phase fluctuations (j> provided that is 
given by 
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1 

v v > 



(18) 



where c 2 = gg is the speed of sound, cf. 0,0- The density 
fluctuations 5g are related to the phase fluctuations <j> via 
5g = -g- l (d/dt + v V)4>. 

If we now calculate the quantum corrections to the 
Bernoulli equation using the effective-action method in 
Eq. by inserting dG^ v jdg or dg^/dg, we obtain 



5 Acs 
do 



;((V0) 2 ) 



(19) 



Clearly, this result differs from the expression (|17|l de- 
rived in the previous Section. Moreover, it turns out 



that Eq. i|17|) contains contributions which are not part 
of the pseudo-energy-momentum tensor (7),„): Even 
though the phonon flux j' in the first line of Eq. I|17|l 

corresponds to (6gV(/) + H.c), which indeed reproduces 
the mixed space-time components of the pseudo-energy- 
momentum tensor (T^„), the phonon density g x contains 
((Sg) 2 ) Icn (where (...) rC n means that the divergent c- 
number \X — XX — ^(0) nas been subtracted already) 



which is not. 



which is part of (T^) rcn but also (0 Z 
(Note that ((j) 2 )rcn cannot be cancelled by the other con- 
tributions.) The expectation value in the second line of 
Eq. l(T7|) is basically ((Sg) 2 ) Icn and thus part of (X^„) re n- 
Finally, the expression ((V%^) ® Vx + H.c.) in the last 
line of Eq. ljT7|l contains (V0 ® V0) rcn which does oc- 
cur in (T Ml/ ) rcn , but also (V5g <8> V5g) ren , which does 
not. One could argue that the latter term ought to be 
neglected in the Thomas-Fermi or local-density approx- 
imation since it is on the same footing as the quantum 
pressure contributions V 2 g and [Vp] 2 , but it turns out 
that this expectation value yields cutoff dependent con- 
tributions of the same order of magnitude as the other 
terms, see Section lYlII below. 



V. CHOICE OF CLASSICAL BACKGROUND 

In order to understand the disagreement between ex- 
pressions (|T7|l and (|T§)) . and to demonstrate the depen- 
dence of the explicit expression for the quantum backre- 
action on the choice of the classical background (see also 
Appendix 0), we will consider two particular alternative 
backgrounds. Firstly, the background is prescribed by 
the density and current generated by the condensed part 
■0c of the field operator ^, cf. Eqs. fTTJ) and (fT2*)) . And 
secondly, we consider (the expectation values of) the full 
density and the velocity potential as a (classical) back- 
ground. 

In the first case, the background current is given by 



Q c v c 



- (t'Vi - H.c.) = j — j 



(20) 



Note that the exact splitting j = j c + j is not unique 
and thus the classical background j c and the fluctua- 
tion j are (in contrast to j) not directly measurable, 
cf. the remark on the Gross-Pitaevskii equation ©. The 
evolution of the phonon flux can be inferred from the 
Bogoliubov-de Gennes equation Q 



d_ . 

di Jx 



I^Qx - 2 V ' (( V ^) ® V X + H.c.) 
~g x V (V + 2g c g) 

-\[{x 2 )v(gm 



H.c. 



(21) 



which enables us to derive the evolution equation for the 
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condensate part via the subtraction j c = j — j 
d 

-Q t 3 c = f(Qc,v c )-2Q c V(gQ x ) 

' "(^) 2 V( 5 <x 2 ))+H.c.] , (22) 

which is the same result as obtained from the modified 
Gross-Pitaevskii equation ||BJ. However, the interpreta- 
tion of the r.h.s. of the equation above as a force density 
is not straightforward, since the equation of continuity is 
not satisfied 



§- t Q C + v-j c = tg^ 2 c w) 2 -(r c fr 



(23) 



I.e., an acceleration of the atoms of the condensate caused 
by a real force could make them leave the condensate or 
induce a change in j c and it is difficult to disentangle 
these two effects. Roughly speaking, choosing the con- 
densate part j c (which is not measurable and not con- 
served) as the classical background would correspond to 
the artificial variable Y of the toy model discussed in Ap- 
pendix [S] whereas the full (measurable and conserved) 
current j corresponds to X. 

In the second case, we use exactly the same variables 
as employed in the effective-action method described ear- 
lier and take the expectation values of the operators cor- 
responding to the density and the velocity potential as 
background 



$ = ($)+0 = b + 0, 
8 = (g) + Sg = g b + Sg . 



(24) 



Note that, in contrast to the full density which is a well- 
defined and measurable quantity, the velocity potential, 
<3E>, is not 14] . It can be introduced via the following 
ansatz for the full field operator 



* = e l 



(25) 



Since $ and g do not commute, other forms such as 
$ = v^e** would not generate a self-adjoint <t (and si- 
multaneously satisfy iJfTij/ = Insertion of the above 
ansatz into the expression for the full current yields 



equals the full density g b 
and deduce 



£>), we remember j = j c +j y 



ghVb = g c v c 



(28) 



i.e., the background velocity is (again up to terms of third 
or higher order) proportional to the condensate velocity. 
Now, considering a condensate which is initially at rest 
v c {t = 0) = and calculating 



Of 



dtj c 



g c d 
g h dt g c + g x 
f(g c = 8b- fto^c = 0) 



g c + g x 



2V( ffft J 



l 

2g~ c 



(C) 2 V( ff (x 2 ))+H.c. 



(29) 



we observe that d t Vb contains the terms g x , (x 2 ), and 

((jc) 2 ), but not ((V0) 2 ). As a result, the quantum cor- 
rections to the Bernoulli equation obtained in the above 
direct manner are in contradiction to the effective-action 
method result SAcs/Sgb- 



VI. FAILURE OF EFFECTIVE-ACTION 
TECHNIQUE 

After having demonstrated the failure of the effective- 
action method for deducing the quantum backreaction, 
let us study the reasons for this failure in more detail. 
The full action governing the dynamics of the fundamen- 
tal fields and * reads 



dt 



(30) 



Linearization \& = i/j c + X yields the effective second- 
order action generating the Bogoliubov-de Gennes equa- 
tions 



* 9 



X 



1 



V' + V + 2g\^ c \ z ) X 



|«)V+H.c. 



(31) 



'g V$ 



(26) 



The quantum corrections to the equation of continuity in 
this formulation can be obtained by inserting the split in 
Eq. I|24|) and neglecting terms of third or higher order 



ghVq 



1 



(SgV(f> + R.c.) = g h V4>b +j" , (27) 



and are in perfect agreement with the effective-action 
method 5A C &/S4>b- 

In order to relate the background velocity V0b = v b 
to known quantities (the background density of course 



Now SA* s /5g c indeed yields the correct backreac- 
tion I122H in the condensate formulation. 

However, if we start with the action in terms of the 
nonfundamental variable $ 



-e[g}-Vg, 



(32) 



with e[g] denoting the internal energy density, the quan- 
tum corrections to the equation of continuity 5A s / 6<fib 
are reproduced correctly but the derived quantum back- 
reaction contribution to the Bernoulli equation SActt/Sgh 
is wrong. 
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Why is the effective-action method working for the fun- 
damental field but not for the nonfundamental vari- 
able $? The quantized fundamental field ^ satisfies the 
equation of motion (3J as derived from the above ac- 
tion and possesses a well-defined linearization via the 
mean- field expansion JSJ. One of the main assump- 
tions of the effective-action method is a similar procedure 
for the variable $, i.e., the existence of a well-defined 
and linearizable full quantum operator $ satisfying the 
Bernoulli equation (for large length scales) 



+ i(Vl>) 2 



h[g] = , 



(33) 



dt 

where h = de/dg + V . The problem is that the com- 
mutator of g and $ at the same position diverges and 
hence the quantum Madelung ansatz in Eq. (|25[) is sin- 
gular. As a result, the above quantum Bernoulli equation 
is not well-defined (in contrast to the equation of conti- 
nuity), i.e., insertion of the quantum Madelung ansatz in 
Eq. into Eq. (QJ generates divergences 

In order to study these divergences by means of a 
simple example, let us consider the generalized Bose- 
Hubbard Hamiltonian 



£ = a H.c.) 



(34) 



where < ij > denote nearest neighbors and hi — ^f^i 
the filling factor of one lattice site i. In the usual contin- 
uum limit, this Hamiltonian generates Eq. (3J where a 
is related to the lattice spacing and the effective mass, 7 
determines g, and V is governed by a, /?, and 7. Insert- 
ing the quantum Madelung ansatz in Eq. <|25|l . however, 
the problem of operator ordering arises and the (for the 
Bernoulli equation) relevant term reads 



1 



H&> — — 



4 ^ 



y/ni(hi + l) (V$)2+H.c. : 



(35) 



with the replacement n, + 1 instead of hi being one ef- 
fect of the noncommutativity. In the superfluid phase 
with large filling n > 1, we obtain the following leading 
correction to the equation of motion 



n lb n 



,(36) 



which depends on microscopic details (such as the filling) . 

By means of this simple example, we already see that 
the various limiting procedures such as the quantiza- 
tion and subsequent mean-field expansion, the variable 
transformation (*&*,*!/) <-> (g, $), and the linearization 
for small fluctuations, as well as continuum limit do not 
commute in general - which explains the failure of the 
effective-action method for deducing the quantum back- 
reaction. The variable transformation ($*,$) <-> (g, $) 
is applicable to the zeroth-order equations of motion for 
the classical background as well as to the first-order dy- 
namics of the linearized fluctuations - but the quantum 
backreaction is a second-order effect, where the afore- 
mentioned difficulties arise Il4l. 



VII. CUTOFF DEPENDENCE 

As mentioned in the Introduction, another critical is- 
sue is the UV divergence of (T M „). Extrapolating the 
low-energy effective action in Eq. to large momenta k, 
the expectation values (Sg 2 ) and (4> 2 ) entering g x would 
diverge. For Bose-Einstein condensates, we may infer the 
deviations from Eq. (IJ at large k from the Bogoliubov- 
de Gennes equations • Assuming a static and homoge- 
neous background (which should be a good approxima- 
tion for large k), a normal-mode expansion yields 



Xk 




Wfc 1 \ A + 



w fe 1 



, (37) 



where d k and dk are the creation and annihilation oper- 
ators of the phonons, respectively, and the frequency LUk 
is determined by the Bogoliubov dispersion relation 
wf, = ggk 2 + fc 4 /4. Using a linear dispersion uj 2 , cx k 2 
instead, expectation values such as (x x) would be UV 
divergent, but the correct dispersion relation implies 

Xk ~ flfcSP/fe 2 + Ofc for large k 2 , and hence (jfix) 1S UV 
finite in three and lower spatial dimensions. Thus the 
healing length £ acts as an effective UV cutoff fc| ut . 

Unfortunately, the quadratic decrease for large k in 
Eq. (22J, Xk ^ o^gg/k 2 + dk, is not sufficient for ren- 
dering the other expectation values (i.e., apart from 
g x and j ) in Eq. (|17|l UV finite in three spatial di- 
mensions. This UV divergence indicates a failure of 
the s-wave pseudo-potential g8 3 (r — r') in Eq. (3} at 
large wave-numbers k and can be eliminated by re- 
placing gS 3 (r — r') by a more appropriate two-particle 
interaction potential Vi nt (r — r'), see 0. Introduc- 
ing another UV cutoff wavenumber fc™* related to the 
breakdown of the s-wave pseudopotential, we obtain 
((Vx f ) Vx + H.c.) - g 2 g 2 fc™* and ( x 2 ) - gg k c s ut . 

In summary, there are two different cutoff wavenum- 
bers: The first one, fc| ut , is associated to the breakdown 
of the effective Lorentz invariance (change of dispersion 
relation from linear to quadratic) and renders some - but 
not all - of the naively divergent expectation values finite. 
The second wavenumber, k^ ut , describes the cutoff for 
all (remaining) UV divergences. In dilute Bose-Einstein 
condensates, these two scales are vastly different by defi- 
nition; because the system is dilute, the inverse range of 
the true potential, which is of order fc^ ut , must be much 
larger than the inverse healing length: 



JL CUt 



jLCUt 

Lorentz ' 



(38) 



Note the opposite relation fc£orentz ^ ^uv i s ver y un " 
natural since every quantum field theory which has the 
usual properties such as locality and Lorentz invariance 
etc., must have UV divergences (e.g., in the two-point 
function) . 

The renormalization of the cutoff-dependent terms is 
different for the two cases: The fc^-contributions can be 
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absorbed by a p -independent renormalization of the cou- 
pling g [l(J, uM > whereas the fc| ut -contributions depend 
on the density in a nontrivial way and thus lead to a 
renormalization of the pressure and the chemical poten- 
tial etc., see the next Section. 



VIII. SIMPLE EXAMPLE 

In order to provide an explicit example for the quan- 
tum backreaction term in Eq. (|17|l . without facing the UV 
problem, let us consider a quasi-one-dimensional (quasi- 
1D) condensate 0,0, where all the involved quantities 
are UV finite. In accordance with general considerations 
|18| . the phonon density g x is infrared (IR) divergent in 
one spatial dimension, therefore inducing finite-size ef- 
fects. Nevertheless, in certain situations, we are able to 
derive a closed local expression for the quantum backre- 
action term Q: Let us assume a completely static con- 
densate v — in effectively one spatial dimension, still al- 
lowing for a spatially varying density g and possibly also 
coupling g. Furthermore, since spatial variations of g and 
g occur on length scales A much larger than the healing 
length, we keep only the leading terms in £/A <C 1, i.e., 
the variation of g and g will be neglected in the calcula- 
tion of the expectation values. In this case, the quantum 
backreaction term Q simplifies considerably and yields 
(in effectively one spatial dimension, where g = gro and 
g = giD now both refer to the ID quantities) 

Q = -V((V X t )Vx)-^V(.g 2 e(2 X t i+tt t ) 2 + ^2^ 



V {ggf' 2 ) + ( ff 5 /V /2 ) + 0(e/X 2 ) 

(39) 



It turns out that the IR divergences of 2(x ! x} an d 
((x ) 2 + X 2 ) cancel each other such that the resulting 
expression is not only UV but also IR finite. Note that 
the sign of Q is positive and hence opposite to the con- 
tribution of the pure phonon density (x x)> which again 
illustrates the importance of the term ((x^) 2 + X 2 )- 

A possible experimental signature of the quantum 
backreaction term Q calculated above, is the change in- 
curred on the static Thomas-Fermi solution of the Euler 
equation lfT3|) for the density distribution (cf. 0, ) 



QlD 



fjt — V VlJ^V 



2tt 



(40) 



with /i denoting the (constant) chemical potential. 
The classical [0(iV)] density profile g c \ — {l^ — V)/giT> 
acquires nontrivial quantum [C(A )] corrections 
Qq = v 1 V" — V /2?r where the small parameter is the 
ratio of the interparticle distance 1/q = 0(1/N) over 
the healing length £ = O(N ). Note that the quantum 
backreaction term gq in the above split g = g c i + gq 
should neither be confused with the phonon density g x 



in g = g c + g x (remember that g x is IR divergent and 
hence contains finite-size effects) nor with the quantum 
pressure contribution oc V 2 y / £) in Eq. (TSJ. 

Evaluating explicitly the change AR of the Thomas- 
Fermi size (half the full length), where fi = V, of a 
quasi-lD Bose-Einstein condensate induced by backre- 
action, from Eq. (|4TJ|l we get AR = — 2~ 5 /' 2 (ujj_/uj z )a s . 
(The quasi-lD coupling constant gin is related to the 3D 
s-wave scattering length a s and the perpendicular har- 
monic trapping lu± by giu — 2a s ui± |l6j.1 In units of the 
classical size R c \ = (3a s A'a;_L/w 2 ) 1/ ' 3 , we have 



1 



1 



AR 



1/3 



2/3 



(41) 



where a z = 1/ describes the longitudinal harmonic 
trapping. In quasi-lD condensates, backreaction thus 
leads to a shrinking of the cloud relative to the classi- 
cal expectation - whereas in three spatial dimensions we 
have the opposite effect 0, ^| . For reasonably realis- 
tic experimental parameters, the effect of backreaction 
should be measurable; for N ~ 10 3 , lu±/uj z ~ 10 3 , and 
a s /a z ~ 1CT 3 , we obtain \AR/R cl \ ~ 1%. 



IX. CONCLUSIONS 

Even though the explicit form of the quantum back- 
reaction terms depends on the definition of the classical 
background, the effective-action method does not yield 
the correct result in the general case (i.e., independent 
of the choice of variable etc.). The knowledge of the 
classical (macroscopic) equation of motion - such as the 
Bernoulli equation - may be sufficient for deriving the 
first-order dynamics of the linearized quantum fluctu- 
ations (phonons), but the quantum backreaction as a 
second-order effect cannot be obtained without further 
knowledge of the microscopic structure (e.g., operator 
ordering). It is tempting to compare these findings to 
gravity, where we also know the classical equations of 
motion only 



1 



Rfiu 7T 9 [iv R — kT[ 



(42) 



which - in analogy to the Bernoulli equation - might 
yield the correct first-order equations of motion for the 
linearized gravitons, but perhaps not their (second-order) 
quantum backreaction. Another potentially interest- 
ing point of comparison is the existence of two differ- 
ent high-energy scales - one associated to the break- 



down of Lorentz invariance fc| ut = fc£"* entz and the other 
'■' (;t to the UV cutoff as well as the question 



LCUt LC 

""UV — A ° 



of whether one of the two (or both) correspond to the 
Planck scale in gravity. 

The dominant 0(£/\) quantum backreaction contribu- 
tions like those in Eq. I|4(J|I depend on the healing length 
as the lower UV cutoff and hence cannot be derived from 
the low-energy effective action in Eq. Q using a covariant 
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(i.e., cutoff independent) regularization scheme, which 
does not take into account details of microscopic physics 
(represented, for example, in the quasiparticle dispersion 
relation). Note that the leading 0(£/A) quantum correc- 
tion to the pressure could be identified with a cosmologi- 
cal term, (T Mi y) = A g^ v in Eq. J2J, provided that A is not 
constant but depends on g and g. (In general relativity, 
the Einstein equations demand A to be constant.) 

As became evident from the remarks after Eq. (|18|l . 
the knowledge of the expectation value of the pseudo- 
energy-momentum tensor (Tn V ) is not sufficient for de- 
termining the quantum backreaction effects in general. 
Even though (T^) is a useful concept for describing 
the phonon kinematics (at low energies), and one may 
identify certain contributions in Eq. (|17|l with terms oc- 
curring in (T Mt/ ) in curved space-times (e.g., next-to- 
leading-order terms in £/A in the gradient expansion), 
we have seen that it does not represent the full dynam- 
ics. Related limitations of the classical pseudo-energy- 
momentum tensor have been discussed in 0]. 

In general, the quantum backreaction corrections to 
the Euler equation in Eq. I|15[) cannot be represented as 
the gradient of some local potential, cf. Eq. (|17|) . Hence 
they may effectively generate vorticity and might serve 
as the seeds for vortex nucleation. Note that this effect 
cannot be observed in the "force" density in the conden- 
sate formulation (|22[1 since j c = p c V0 c holds and thus 
introducing vorticity into the condensate requires creat- 
ing condensate holes with g c — (vortices), which is be- 
yond the regime of applicability of our linearized analysis. 
Therefore, vorticity (in the full-j formulation) can only 
be generated by extracting atoms from the condensate. 

In contrast to the three-dimensional case (see, e.g., 
[TTL IT^I. the quantum backreaction corrections given 
by Eq. i|39[) diminish the pressure in condensates that 
can be described by Eq. (0J in one spatial dimension 
(quasi-lD case). This is a direct consequence of the so- 
called "anomalous" term ((x^) 2 + X 2 ) m Eq. which 
- together with the cancellation of the IR divergence - 
clearly demonstrates that it cannot be neglected in gen- 
eral. We emphasize that even though Eqs. (|39[) - l|41[l de- 
scribe the static quantum backreaction corrections to the 
ground state, which can be calculated by an alternative 
method 0] as well, the expression in Eq. I|17l) is valid 
for more general dynamical situations, such as expand- 
ing condensates. The static quantum backreaction cor- 
rections to the ground state can be absorbed by a re- 
definition of the chemical potential fj,(g) determining a 
(barotropic) equation of state p(g)', this is however not 
possible for the other terms in Eq. I|17l) , such as the quan- 
tum friction-type terms j®v etc., occurring in more gen- 
eral dynamical situations. 

We have derived, from the microscopic physics of dilute 
Bose-Einstein condensates, the backreaction of quantum 
fluctuations onto the motion of the full fluid and found 
possible experimentally observable consequences. We ob- 
served a failure of the effective-action technique, Eq. 



and a cutoff dependence of the backreaction term due to 
the breakdown of covariance at high energies. Whether 
similar problems beset "real" (quantum) gravity remains 
an interesting open question. 
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APPENDIX A: TOY MODEL 

To give an illustration of the influence of the choice 
of coordinates describing the classical background, let 
us consider a very simple quantum system, the one- 
dimensional harmonic oscillator with the frequency ft. 
The equation of motion for the position operator X reads 

^X + tfX = 0. (Al) 

Of course (as is well known from the theory of coher- 
ent states), splitting up this position operator X into a 
classical (background) part X^ = (X) and a quantum 
fluctuation part SX with (5X) = via X = X\, + SX 
yields the same equation of motion for both parts sepa- 
rately, i.e., there is no backreaction at all. 

However, if we artificially introduce another variable 
Y via X — Y 2 /2, its (classical) equation of motion reads 

Now the same procedure Y = + SY would yield a non- 
vanishing quantum backreaction term. Therefore, one 
must be careful when comparing different expressions 
for the quantum backreaction, since the explicit form 
may depend on the choice of splitting the system into 
a classical background and quantum fluctuations (e.g., 
X = X h + SX versus Y = Y b + SY). This illustrates the 
importance of working with measurable quantities such 
as X or j and g. Exploiting the simple example a bit 
further, an an-harmonic oscillator d 2 X + fl 2 X = gX 2 
yields 

J^(A) + Q 2 (X) g(X) 2 = g(SX 2 ) , (A3) 
i.e., a real quantum backreaction term. 
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